New Exact Solutions of a Generalized Shallow Water Wave Equation 



o 

(N 



Bijan Bagchi 1 * Supratim Das 1 Wd Asish Ganguly 2 * 



department of Applied Mathematics, University of Calcutta, 
92 Acharya Prafulla Chandra Road, Kolkata-700009, India 
2 Department of Mathematics, Indian Institute of Technology, Kharagpur 721302, India 

February 16, 2010 



GO . 

£H \ Abstract 

In this work an extended elliptic function method is proposed and applied to the generalized shallow water 
i_jwave equation. We systematically investigate to classify new exact travelling wave solutions expressible in 
^ terms of quasi-periodic elliptic integral function and doubly-periodic Jacobian elliptic functions. The derived 
j> new solutions include rational, periodic, singular and solitary wave solutions. An interesting comparison with 
CN the canonical procedure is provided. In some cases the obtained elliptic solution has singularity at certain 
^ region in the whole space. For such solutions we have computed the effective region where the obtained 
C<~) solution is free from such a singularity. 

(N : 

o : 

ey words: Shallow Water Wave Equation, Integrable Systems, Travelling Waves, 
Jacobian Elliptic Function, Rational Solution, Singular Solution. 

•i-H , 

X 

PAfeS number(s): 02.30.Jr, 02.30.1k, 05.45,-a 



* bbagchi 1 23 @rediffmail.com 
t supratimiitkgp @ gmail.com 

taganguly@maths.iitkgp.emet. in, gangulyasish@rediffmail.com 



1 



1 Introduction 

Seeking analytical as well as numerical solutions of nonlin- 
ear systems has continued to attract attention through the 
last few decades HI El |3l @J 0. Mathematically, nonlin- 
ear equations do not normally have solutions which would 
superpose making the systems they represent rather com- 
plicated and difficult to analyze OQ. On the other hand, 
an extensive study of a number of nonlinear systems has re- 
vealed that there do exist solutions which are not only inter- 
esting in their own right but have a wide range of applica- 
bility [8l[9]|. To generate exact solutions and to understand 
their properties several important techniques have been de- 
veloped such as the inverse scattering approach IfTOll . Lax 
pair formulation [11 J, Backlund transformations ][T2l and 
Hirota's bilinear method IPT31 [Pfll . Nonlinear equations 
have also been shown to arise from a boundary value prob- 
lem [15] possessing hierarchy of conservation laws [fT6ll . 

An important class of solutions of nonlinear evolu- 
tion equations is concerned with those of the travelling 
waves that reduce the guiding partial differential equation 
of two variables namely, x and t to an ordinary differ- 
ential equation of one independent variable z = x — ct 
where c(g K — {0}) is a parameter signifying the speed 
with which the wave travels either to the right or left. A 
number of methods has been employed in the literature to 
obtain the travelling wave solutions of various types. Of 
these the tanh-method stands out as one of the very ef- 
fective tools for solving certain classes of nonlinear evo- 
lution and wave equation |fl7ll . Others include the varia- 
tional iteration method for a class of linear and nonlinear 
equations IfTSTl . homogeneous balance method |[T9l , the hy- 
perbolic method [20], the trigonometric method [21], Dar- 
boux's transformation [1221 . F-function method (231, G'/G- 
expansion method [|24l . a unified algebraic method [|25l 
and more. 

Recently the existence of travelling wave solutions for 
the generalized shallow water wave (GSWW) equation 

m 

du 

u xxxt + au x u xt + (3u t u xx -u xt -u xx = 0, (ui = — ) (1) 

with a, (3 E R — {0} has been noted. The derivation of 
(OQ) follows from the classical water wave theory with the 
aid of Boussinesq approximation. In an interesting review, 
Clarkson and Mansfield [|27l considered the various clas- 
sical and non-classical reductions of the GSWW equation 
wherein they also investigated the Painleve tests to examine 
the complete integrability of (OQ) which holds if and only if 
a = f3 or a = 2/3. They explored some simple and non- 
trivial family of solutions of (OQ) while in ll28l a class of 
exact travelling wave solutions were obtained by making 
use of the homogeneous balance method and a modified 



hyperbolic method. 

In this article we plan to study the travelling wave solu- 
tions 

u(x, t) = u(z) , z = x — ct , (2) 

of the GSWW equation (OQ) in a more general framework. 
We show that the canonical procedure namely, the classi- 
cal technique [1291 , based on an integration within a suit- 
able range not only recovers both known periodic (includ- 
ing elliptic and trigonometric) and non-periodic hyperbolic 
solutions, but also unfolds several new elliptic and ratio- 
nal solutions. However such a constructive method has its 
limitations in the sense that the range of validity of the ob- 
tained solutions is fixed by the analysis of the zeros of the 
governing cubic polynomial in u'(z) beforehand and so, in 
practical applications, a particular solution proves difficult 
to implement. 

One of our aims is to go beyond this standard path by 
proposing an extended method based on the use of Jaco- 
bian elliptic functions. We refer to it as an extended elliptic 
function (EEF) method. In this context the works of Refs 
[[30l and OTl worth to be mentioned which used the ellip- 
tic functions for a class of nonlinear evolution equations. 
These works consider a finite series expansion in the form 

rn 

v(z)=u'(z) = Y,^F l (z), (3) 

i=0 

to generate travelling wave solutions. In ([3]) the positive in- 
teger m is determined by balancing the coefficients Oj's af- 
ter substituting the expansion in the given nonlinear equa- 
tion and comparing the relevant powers. Our EEF method 
generalizes this approach by including the negative powers 
of F thereby addressing the full-range series as given in the 
following 

m 

v(z) = a i F % z ) • ( 4 ) 

i=—m 

In Equation ©, F is an unknown functional constrained to 
satisfy the relation 

( F i)i = (l + £lF 2)( £2 + £3F i), (5 ) 

where the prime denotes the derivative with respect to z 
throughout the text. The important feature of the EEF 
method is that it includes not only the positive integral 
powers of F but also the negative ones. Our motivation 
for keeping lower negative degrees of F comes from the 
fact that in typical solvable systems like the Korteweg de 
Vries the related spectral problem allows similar exten- 
sions [[321 [331 . It has been shown that these lower equations 
may arise 1(331 as a necessary part of an extended symme- 
try structure [34] while it is known [[331] that the underlying 
commutator generates the whole set of symmetry operators 
for the system. 
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It is not difficult to see that the biquadratic integral 061 
emerging from the constraint © facilitates expressing the 
travelling wave solutions in terms of the Jacobian elliptic 
functions pq(z — z ,k), where the symbol pq represents 
twelve such distinct functions sn,cn,dn,ns(= l/sn),cs(= 
cn/sn) etc. The novelty of the EEF method is that it 
opens up a broad spectra of new travelling wave solu- 
tions that include the already known ones (see for exam- 
ple, J22 ED [371 EH ED). We shall explicitly provide the 
solutions generated from the EEF method for some special 
selections of the integral parameters. 

This article is organized as follows. In Section [21 we 
derive the reduced ordinary differential equation for the 
GSWW equation and then go for the classical technique 
to obtain exact analytical solutions. In Section [3] we turn 
to the constraint © to take up the construction of our new 
procedure and then categorise different types of solutions 
being generated from the EEF method by appealing to par- 
ticular choices of the integral parameters. It will be seen 
that the solutions in the appropriate limits reduce to pre- 
viously known results. Finally, in Section |4] we present a 
summary of our results. 

2 Canonical procedure for the 
GSWW equation 

To start with, we substitute the form © into the GSWW 
equation © that results in a fourth order ordinary differen- 
tial equation 



cu"" + (a + /3)cu'u" + (1 - c)u" = . 



(6) 



Let us note that the parameters a, (3 controls the strength 
of non-linearity of above equation. Without loss of gen- 
erality, one can leave one of them, say a, arbitrary while 
the other parameter (3 may be varied to tune the strength of 
non-linearity appropriate for the concerned physical situa- 
tions. Three particular values of (3 namely j3 = ±a, a/2 
are known to correspond completely integrable system. In 
this article we shall explore the solutions of Equation © 
leaving both of them arbitrary except for /3 = —a as it 
points to the linear system and hence is well-known. It is 
straightforward to obtain a second integral of © given by 

3cM //2 + c(a + /3y 3 + 3(l-cy 2 + 6cciu' + 3cc 2 = 0, (7) 

which may be expressed as a first order equation in terms 
of the variable v(z) = u'(z): 



gration are denoted by the symbol Cj and will not be explic- 
itly mentioned further. In © the monic cubic polynomial 
"P^v) is given by 



m / \ i 3(1 — c) 2 6ci 3c 2 
c(a + p ) a + p a + p 



(9) 



For simplicity let us assume that all the three roots of 7 3 (v) 
are real and focus on the arrangement V\ > v 2 > v 3 . 
We consider first the non-degenerate case of distinct roots 
for a + (3 < 0. A formal integration of © in the range 
[(vj,Zj),(v,z)], where the point Zj is to be determined from 
the transcendental equation v(zj) = Vj, yields 



dz = ± 



dw 



n;=i(^-^-)] i/2 



(10) 



Now care is to be taken to choose the root Vj,j = 1, 2 or 3 
that defines the interval of validity of the solution. It is easy 
to see that the reality condition for v(z) gives rise to two 
different cases namely either v 3 < v < v 2 or v 1 < v < oo. 
We address them by turn below. 



Case 1: 



v 3 < v < v 2 (a + (3 < 0) 



Keeping in mind that the range of integration is between 
[t> 3 , v] and [za, z], the following change of variables w — >■ (p 
proves useful 



w 



V3 + (v 2 - v 3 ) sin Lp . 



(11) 



The standard definition for Jacobian elliptic function then 
expresses v(z) in a closed analytic form 

v(z) = u\z) =v 3 + (v 2 - v 3 )sn 2 (u(z - z 3 ), k) , (12) 



where v 



y/[\a + (3\{ Vl - i> 3 )/3]/2 and k 2 



{v 2 



v 3)/( v i ~ u 3)(0 < k 2 < 1) is the modulus of the ellip- 
tic function BOl . In terms of the elliptic integral function 
of the second kind E(ipi,k) through the relation siny?£ = 
sn(u (z — zi),k), defined by the integral 

E(tp e ,k)= VI - k 2 sin 2 a da, £=1,3,(13) 
Jo 

we obtain u(z), after integrating (fT2l) . in the final form 



u(z) = u 3 + 



v 2 -k' 2 v 3 v 2 -v 3 

\z-z 3 ) —E(ip 3 ,k). (14) 



A- 2 



uk 2 



,'2 



a + (3 



a + (3 



n< 



V — Vi 



(8) 



Note that in Equation © and also elsewhere in the text, the 
arbitrary constants appearing through the process of inte- 



rn the above solution and also in the following Uj is defined 
in a natural way as Uj = u(zj),j = 1,2,3 while k' 2 = 
1 — k 2 is the complementary modulus of elliptic functions. 

It may be pointed out that (fT4l) is a new travelling wave 
solution. Two interesting limits namely k — > 1 and k — > 
may be of interest to notice that lead to the existence of a 
double root of 7 3 (v). For k — > l-o, we must let v\ — > v 2 +o, 
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but then v 3 cannot be allowed to coincide with v 2 as this 
will imply that the limiting solution has no range of va- 
lidity. Along the same reasoning we conclude that the 
limit k — > is disallowed for this case. Noting that 
E(ip 3 , k — > 1) = tanh[z/(;z — z 3 )} the solution (fl4b thus 
reduces to the hyperbolic form for k — > l-o 



Case 4: 

v{z) 
u(z) 



-oo < v < v 3 (a + (3 > 0) 

ui - (t>i - v 3 )dc 2 {v{z - z 3 ), k') , 

U 3 + V 3 { Z - Z 3 ) + [£(<^, fc') 



(23) 



u(z) = u 3 + v 2 (z- z 3 ) 



v 2 - v 3 



V 



tanh[^(z — z 3 )\ . (15) 



-sn(i/(z - z 3 ), k')dc(v{z - z 3 ), k')} , 
which comes from the following substitution in ([10 



Case 2: Vj < v < oo (a + /3 < 0) 

Clearly the range of integration now lies between [v\ , v] 
and [zi, z\. Employing the substitution 



w 



Vi — (v i — w 3 )cosec 2 <^ . 



(24) 



f3 + (vi — v 3 )cosec 2 Lp 



(16) 



we obtain the following expression for v(z) 

v{z) = u'(z) = v 3 + (vi - v 3 )dc 2 (u(z - zi),k) . (17) 

The above equation, on further integration, yields another 
new travelling wave solution 



u{z) 



Mi + vi(z- zi) - — — — [E(tpi, k) 



It is interesting to note that in contrast to the Cases 1 and 2 
for a + (3 < 0, the known hyperbolic and periodic solutions 
are recovered for a + (3 > (Cases 3,4) in the complemen- 
tary limit of the elliptic solutions (12TT) and (|23l . 

Finally let us deal with the degenerate case when 7 3 (v) 
has a triple zero, i.e. v i = v 2 = v 3 = v . From the cubic © 
it is clear that we have to choose c 2 = — 2u ci/3 giving the 
multiple root v = (c — l)/[c(a + /?)]. The range of inte- 
gration will be — oo < v < v or v < v < oo according as 
a + f3 ^ 0. The singular rational solution can be expressed 
in a compact form as follows 



v 



-sn(u(z-z 1 ),k)dc(iy(z-z 1 ),k)} . (18) ? 3 (v) = (v-v ) 



3 . 
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In above we have used the relation snf^ + i^, k) = cd(z, k), 

where K(k) = J^ 2 da/ y/l — k 2 sin 2 a is the quarter real 
period of elliptic sine function. Let us now take the limit 
v 2 — > v i-o keeping v 3 < v 2 that makes the travelling wave 
solution linear in spatial and time coordinates 



k — > 1-0 : u(x, t) = {u\ — V\Z\) + V\(x — ct) 



(19) 



On the other hand, letting v 3 — > v 2 -o where v 2 < V\, we 
recover the periodic solution 

k — > : u{z) = U\ + v 2 (z — Z\) H — tan[z/(2; — z\)\ . 

(20) 

The case for a + f3 > will give the result in comple- 
mentary frame which can be derived along the same way 
as before. Without giving the details of the procedure, we 
are providing below the final solutions. 



Case 3: 



v 2 < v < vi (a + (3 > 0) 



v(z) = Vx - (vi - v 2 )sn 2 (v(z - Zt), k') , 

( \ i ^2 - k 2 v i Vl-^2 , , 

=ui + — (z - Zl ) + — E((p x , k ) ,(21) 



>k' 2 



where sin<£>- = sn(v(z — zi),k'),i = 1,3. To derive the 
above we made the following substitution in (TTOl) 



u = v {z-c 3 )^ — hc 4 , 

a + p z — c 3 

(25) 

which was obtained in 11281 from a different approach. 

We have thus obtained several new solutions by apply- 
ing the classical technique based on an analysis of the zeros 
of a cubic polynomial. Interestingly in the appropriate lim- 
its all the known solutions are retrieved. It is clear that for 
practical applications of the solutions u(z) derived in this 
section, a knowledge about the position of zeros of u"(z) is 
required which somewhat weakens their utility. In the next 
section we propose a new method that is free from such a 
limitation. 



3 The EEF method for the GSWW equation 

Here instead of focussing on the zeros of u"(z) which are 
unknown we fix the zeros of u'(z). Since u" has three zeros 
[see Equation ([8]) of Sec. |2), it is obvious that u' has four 
zeros. This motivates us to propose the following construc- 
tion. 

3.1 Construction 

Let us turn to the fourth-order ordinary differential equa- 
tion © which we integrate to write 



2cv" + c(a + /3)v 2 + 2(1 - c)v - c 5 = 0. 



(26) 



w 



Vi + {v 2 - Vi) sin ip . 



(22) 
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Table 1 : The solutions of the non-linear system d29j > for the expansion 
parameters aj , j = 0, 1 , 2 are given. In the second column, uj = E\Ei + 

£3- 



Table 2: The integration constant C5 and a_i, a_2 satisfying d29l > are 
provided where (Aj - w 2 )/(£i£ 2 £ 3 ) = 12/(— 4)- 7 ' -1 , for j = 1,2 and 
yi 3 =ca (a + /3) + 2(l-c). 



Class 


(a + /3)a 


ai 


(a + /3)a 2 


Class 


a_i 


(a + /3)a_ 2 


c(a + (3)c 5 


1 


1 -c- 1 -4w 





-12ei£ 3 


1 





-12e 2 


-(l-c) 2 + 16c^i 


2 


1 -c" 1 - 4w 








2 





-12e 2 


-(l-c) 2 + 16c 2 il2 


3 


1 -c- 1 - 4w 





-12ei£ 3 


3 








-(l-c) 2 + 16c 2 il2 


4 


Arbitrary 
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ca A 3 (a + P) 



We look for a formal solution in the form 

m 

v(z)=Y, a * Fl i ( 27 ) 

i=— m 

where the generating functional F(z(x,t)) satisfies the 
constraint ©. Note that the inclusion of negative powers 
of F, in general, creates pole of u' triggering the presence 
of singular solutions in some situations. However, we ex- 
plicitly show that such singular solutions are physically ac- 
ceptable in a restricted domain of space labeled by (x, t,u). 
Substituting the expansion (1271) into Equation (1261) points to 
m = 2 and so the Lorentz-like expansion of v(z) reads 

v(z) = a + a 1 F + a 2 F 2 + ^ + ^. (28) 

We next compute the expansion parameters that 
needs term-by-term balancing of the coefficients of each 
powers of F to zero. Somewhat involved but straightfor- 
ward algebra leads to the following relations 

F ±4 : a ±2 [12A±+(a + /3)a ±2 ] =0, [A+ = e 1 e 3 ] , 
F ±3 : a ±1 [2A±+(a + /3)a ±2 ] = 0, [A_ = e 2 ] , 
F ±2 : 8c(£i£: 2 + e 3 )a± 2 + (a + /3)c(a| 1 + 2a a ±2 ) 

+ 2(1 - c)a± 2 = , (29) 
F ±x : c(ei£: 2 + e 3 )a ±1 + (a + /3)c(a a ± i + a Tl a± 2 ) 

+ (1 -c)a±i = 0, 

F° : 4c(£ 2 a 2 + £i£ 3 a_ 2 ) + (a + /3)c(a^ 
2 

+ 2 ctia_i) + 2(1 — c)a + c 5 = . 

By exploiting (|29l) we can derive four classes of solu- 
tions for aj and c 5 which are summarized in TableQ]and Ta- 
ble [2l From a previous work we already know that twelve 
different choices exist (see Table I in lT36lO for the zeros 
of F in right-hand side of © that leads to different repre- 
sentation for F(z) in terms of Jacobian elliptic functions. 
Further for each of them one gets four classes of travelling 
wave solutions of GSWW equation (QQ) corresponding to 
the solutions furnished in Table [TJ and Table [21 Our new 
method is a generalized procedure and encompasses the 
previously known solutions as special cases. In the next 



subsection we provide, as illustrative cases, four classes of 
solutions for three particular choices of integral parameters 
in Equation ©. 

3.2 Class 1-4 solutions for particular selections 
of integral parameters 

To get explicit forms of Class 1-IV solutions obtained in 
EEF method, it remains to choose the integral parameters 
£j,j = 1,2,3 in the constraint F' 2 = (1 + e 1 F 2 )(e 2 + 
e 3 F 2 ). At first let us consider some degenerate selections 
(i. e. taking a pair of double zeros of F') leading to hyper- 
bolic, trigonometric and linear solutions. 

• Algebro-hyperbolic Waves 

The choice (£i,£ 2 , £3) = (—1,1,-1) gives F(z) = 
tanh z that generate the following solutions modulo a 
constant 

u = (do + a 2 + ci- 2 )z — a 2 tanh z — a_ 2 coth z . (30) 

The above solution reduces to linear form for Class 4 
while singular term disappears in Class 3 solutions. 

• Algebro-trigonometric Waves 

Choosing (ei ) e 2) e 3 ) = (1,1,1), we get F = tan z 
which gives following solutions apart from an inessen- 
tial factor 

u = (do — a 2 — a_ 2 )z + a 2 tan z — a_ 2 cot z . (31) 

Note that such types of algebro-hyperbolic and algebro- 
trigonometric solutions are new and of interest. By choos- 
ing suitable value of wave speed c the algebraic term can 
be removed causing the reduction to known forms obtained 
in [|28l . Now using a canonical procedure we have already 
obtained [ see (fT5T) & (|20l) of Sec. O such types of solu- 
tions in the extreme limits of modulus parameter of elliptic 
functions. This means that more general forms of solutions 
can be generated from EEF method provided F' has four 
simple zeros. 

We already mentioned that twelve different selections 
are possible for the simple zeros of F' in the constraint © 
leading to closed analytic expressions for F(z) in terms of 
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doubly periodic Jacobian elliptic functions. Here we pro- 
vide explicit forms of the wave-solutions for the three rep- 
resentative selections, two of which produce singular wave- 
solutions of Class 1 and 2. Note that given non-linear pa- 
rameters a, (3 we fix the the integration constant c 5 of Equa- 
tion (|26)) for a particular value of wave- speed c through a 
relation as dictated by the last column of Tableland hence 
that will not appear in the final expressions of the solutions. 




Figure 1: Class 1 travelling wave for Selection I. The primary inputs: a = 
0.5, /3 — 0.1, k — 0.7 and the scaling parameters are uo — zo — 0. The 
humps are the signature of singularity. 



Selection I 

Let us choose the zeros of F' as ±l,±l/fc(0 < 
k 2 < 1) which correspond to the selection of triplet 



(£i, £2, £3) - 
z such that F(z ) 

of OH) 



-1, 1, -k 2 ) or (-k 2 , 1, -1). Choosing 



0, one then obtains (see Table I 

F{z) = sn(z - z ) . (32) 
Thus we are led to four classes of solutions of which 



the first two (j 
lutions : 



1, 2) correspond to the singular so- 



Class j : 



,(i) 



a +a 2 sn (z-z )+a_ 2 ns (z—Zq), 

(33) 

where c^s are obtained from 1st and 2nd rows of Ta- 
ble Q] & Table [2l The corresponding travelling wave 
solutions are obtained from the integration of (|33l) 
and involves a singularity at z = z as is evident from 



the formula [40| 



Z-ZQ 



ns 2 r 



4 I dr 



z- z + 
z )ds(z - 



1 



z - z 
z )- (34) 



-E(& ,k) -cn(z 

In the expression (|34l) and also in the following di is 
defined by sin^ = sn(z — zg, k) for I = 0, 1. In this 
context it is worth mentioning that in the limit k — > 
1, the travelling wave reduces to coth-type solution 
which also has a singularity at z = z . This solution 
was derived from a different approach in Ref. |[28Tl . It 
is well known that such a singular solution may serve 
as a model for physical phemenon of so-called "hot- 
spots" eh nana. 

In the following our aim would be to find regions 
where the travelling wave solutions (j = 1,2) cor- 
responding to (1331) are regular. We use the follow- 
ing infinite trigonometric series expansion [|40l of 

ns(z — Zo, k) = l/sn(z — z , k) : 



ns z 



zo,k) 

2tt 
~ ~K 



n=0 



7T 

2K 

g 2n 



cosecr 



q 



2n+l 



sin[(2n + l)r] ,(35) 



where q = exp[— irK'/K] is known as nome of the 
elliptic functions and r = n(z — z )/2K, K'{k) = 
K(k'), the quantities K{k) and k have already been 
defined in Sec. [2l Denoting the singular term in the 
travelling wave solutions by the symbol So, we see 
that as z — > z , the asymptotic expression of its 1st 
term will be 



S Q (z -> z ) 



71 

2K° 



rcosecr 



(36) 



We will now use the well known expansion for cosecr 
convergent in the region \r\ < n 



cosecr 



-l) n - 1 2(2 



2n-l 



^)B 2n 



(2n)\ 



-T 



.2n-l 



(37) 

where {£>„} are Bernoulli numbers whose first few 
members are B = l,B t = —1/2, B 2 = 1/6, B A = 
-1/30 etc. 

Combining (1361) and (1371) . it is straightforward to con- 
clude that near z ~ z , the travelling wave solution 
behaves like (z — zo) 1 i- e - the point z = zq is a sim- 
ple pole. Note that the region of convergence of the 
series (1371) for cosecr is 



|r| < 7T \z - z \ < 2K (38) 
The above equation speaks of two parallel planes 

x - ct - z ± 2K = , (39) 
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which are normal to wave surface. These two planes 
are shown in Fig.[T][see two extreme strips orthogonal 
to surface u = u(x, t)]. It is thus quite expected that 
wave surface will be irregular along these two planes. 
This fact can be appreciated from the humps along 
the strips in Fig. [TJ showing sudden rise of the value 
u = u(x — ct). Further there must be a singularity at 
z = zq. From elementary analysis we know that in 
t-x plane the point z = x — ct approaches point z 
along any curve. Hence there will be a circular disc in 
the plane z = z centred at z and of arbitrary small 
radius where solution will show irregularity. This is 
reflected by a single localized hump in Fig. [T] on the 
plane lying parallel and midway between two extreme 
planes. 

To use such a solution we thus need an exact knowl- 
edge about the position of parallel humps. These 
may be exactly computed by the solid angles 6 X = 
cos -1 |1/a/1 + c 2 | and 6^ = cos -1 \c/y/l + c 2 |. These 
two angles in turn depend on the parameter c which 
determines the wave speed. Hence very fast and very 
slow waves have different regions of humps. The core 
singularity at z = z will however be present in same 
region for all types of waves which may be removed 
by making a hole of arbitrary small radius on the plane 
z = zq around point z . This concludes our detailed 
analysis about singularity of Class 1 and 2 solutions 
for Selection I which hopefully provides a sufficient 
basis of the solution regarding its potential use. 

The other two classes correspond to the non- singular 
solutions 

Class j : v {j) = a + a 2 F 2 , j = 3, 4 . (40) 

The resulting travelling wave solutions are of the 
same form as the solutions (PT41) and (TT9b obtained in 
Sec 12 

w (i) = Uo + ( fl0 + ^)(z~ zo) - ^E(# , k) , (41) 

where a^s are to be computed from last two rows of 
Table [TJ and Table [2l Note that Class 4 solution is 
linear in x and t, since a±\ = a± 2 = (see the last 
rows of the tables). 

• Selection II 

Let us now choose pairs of purely real and of purely 
imaginary roots of F' respectively as ±l/k'(0 < 
k' 2 = 1 — k 2 < 1) and ±i/k which come from the 
selection (e 1 ,e 2 ,e 3 ) = (—k' 2 ,l,k 2 ) or (k 2 , 1, — k' 2 ). 
One then obtains the following representation of F(z) : 

F(z) = sd(z-z ,k), (42) 




Figure 2: Class 1 travelling wave for Selection III for the same primary inputs 
and scalings are ui = z% = 0. The smoothness of the wave-surface shows the 
solution is non-singular. 

leading to a new solution of u(z) . The explicit expres- 
sions for Class 1 and Class 2 are given by 

u {]) = u + Qj(z - z ) - (3 - j)a_ 2 -E($o, k) 
+ (2-])k 2 a- 2 sn(z-zo, ^c^z-zq, k) + S (43) 

for j = 1,2, where we have abbreviated qj as Qj = 
a + (3 — j)k' 2 a^2- The non-singular solutions corre- 
sponding to Class 3 and 4 solutions read 

uU) = u o+ (ao zo) + ^E^o, k) 

- ^|sn(2; - z , k)cd(z - z , k) , j = 3, 4 . (44) 

The main feature of Selection II is that the solutions 
coming from that are nearly periodic in the whole do- 
main due to the additional sn and cd-term. In con- 
trast the solutions generated from Selection I are only 
quasi-periodic, the periodic behaviour is observed for 
singular solutions which are prominent near z ~ z$. 

• Selection III 

We just saw that the Class 1 and 2 solutions generate 
singular solutions owing to the generating functional 
F(z(x,t)) possessing either a zero or a pole in the 
region of validity. It is indeed possible to generate 
non-singular solutions for these classes, which share 
similar qualitative behaviour with the previous coun- 
terparts, if F can be chosen to have neither zero nor 
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a pole in the finite part. Below we provide such an 
example, which are the only one among the set of el- 
liptic functions. 

We choose the zeros of F' as ±1, ±l/k' correspond- 
ing to the selection (e^, e 2 , £3) = (—k' 2 , — 1, 1). Then 
we have the following representation of F(z) (see Ta- 
ble I of 11361) 



References 



F(z) = nd(z — Zi, k) 



(45) 



where z\ is fixed as F(^) = 1. It is now trivial to 
check that the generated solutions for each of the four 
classes remain non-singular 
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U) 
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Ci) 



- ui + a (z - z\) + (3 - j)a_ 2 E(^i, k) 
(2— j)k 2 a_ 2 sn(z — Zi, k)cd(z—Zi, k) ; (46) 
= Mi + a (z - zi) + a 2 E(di, k)/k' 2 , 



k 2 a 2 sn(z - zi, k)cd(z - z 1} k)/k' 2 , (47) 



where j = 1,2 for Equation (1461) and j = 3,4 for 
(|47|) . Fig. |2] depicts a Class 1 wave guided by (|46l) . 

Other possible selections from Table I of 11361 will gen- 
erate many such new elliptic travelling wave solutions of 
GSWW equation. Let us mention that the degenerate selec- 
tions leading to hyperbolic and trigonometric type waves 
can be obtained from elliptic solutions in k — > 0, 1 limits. 

4 Summary 

In this article we proposed an extended method to generate 
a rich class of doubly-periodic elliptic travelling wave solu- 
tions of the GSWW equation. A systematic classification is 
given for the solutions to exhaustively utilize the strength 
of the proposed method. We also discussed a canonical 
procedure for generating travelling wave solutions. The 
problem with such solutions is that these require the pre- 
knowledge of zeros of second derivative of the solutions. 
The proposed EEF method removes this difficulty by fix- 
ing the zeros of the first derivative of the solutions and 
leads to a wide range of travelling waves which include 
the kink type solitary wave lfl6l . sinusoidal type periodic 
solution ll27l and a rational solution (" 112910 . Class 1 to 4 
classify the various types of solutions whose explicit forms 
are noted for three representative selections of integral pa- 
rameters defining the generating functional. Classes 1 and 
2 produce singular solutions which has important applica- 
tion in modelling the physical curcumstances of formation 
of hot-spots iflTl l42l l43l . We have exactly computed the 
positions of such blow-up of solutions and graphically il- 
lustrated in the Fig. [TJ To the best of our knowledge such 
an analysis about the region of blow-up for singular elliptic 
solutions has not been done previously in the literature. 
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